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SECTION  I 
INTRODUCTION 

The  requirements  for  mass  spectrometric  analysis  of  gas  samples 
are  well  known  and  have  been  described  in  detail  (Reference  1).  Analysis 
of  glow  discharges  in  molecular  gases  or  gas  mixtures  is  more  difficult 
because  of  composition  change  due  to  electrochemical  activity.  One 
example  is  that  of  the  sealed-off  carbon  dioxide  laser  which  is  the 
subject  of  this  report.  Mass  spectrometric  measurement  of  steady-state 
component  concentration  in  the  CO2  laser  discharge  has  been  shown  to  be 
feasible  (Reference  2).  Measurement  of  reaction  rates  and  transient 
concentration  is  also  possible  (Reference  3).  The  evidence  indicates 
that  valuable  information  concerning  the  plasma  chemistry  of  the  discharge 
can  be  obtained  with  the  aid  of  a  mass  spectrometer  (References  3-22). 

The  dissociation  of  C02  is  a  major  problem  in  sealed-off  C02  lasers. 
The  carbon  monoxide  formed  competes  with  the  CO^  for  vibrational  excita¬ 
tion;  and  the  atomic  and  molecular  oxygen  formed  affects  the  stability 
of  the  discharge  (Reference  14).  Dissociative  attachment; 

C02  +  e"  -*■  CO  +  O'  (1 ) 


and  collisional  dissociation: 


C02  +  e  ■*  CO  +  0  +  e  (2) 

both  contribute  to  the  dissociation  process  (References  5  and  8).  The 
dissociation  of  C02  seems  to  obey  the  general  reaction 

C02  ■>  CO  +  1/2  02  (3) 

with  an  equilibrium  constant  that  varies  with  the  diameter  of  the  tube 
(Reference  9)  and  at  pressures  less  than  0.5  torr  (Reference  6).  The 
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a 

dissociation  coefficient,  which  is  defined  as  —  where  a  is  the  number  of 

p 

dissociations  occuring  per  electron  per  cm  of  path  length  and  p  is  the 
gas  pressure,  also  increases  with  reduced  field  strength  (E/p)  (Reference 
7). 


CO2  dissociation  can  be  reduced  by  the  addition  of  to  the  dis¬ 
charge  (References  10  and  11).  The  acts  as  a  catalyst  to  reform  the 
CO2,  and  the  equilibrium  constant  for  reaction  (Equation  3)  is  found  to 
be  linearly  dependent  on  the  ^  concentration  (Reference  5). 

The  major  neutral  contaminants  in  a  (^/^/He  discharge  appears  to 
be  CO,  02,  NO,  N^O ,  and  N02  (References  3  and  4).  CO  and  0?  concentrations 
are  comparable  to  the  concentrations  of  the  initial  fill  gases  ^/No/He, 
and  the  nitrogen  oxides  are  in  concentrations  less  than  1%.  Although 
the  concentrations  of  the  nitrogen  oxides  are  small ,  the  laser  gain  is 
affected  by  them  (Reference  12). 

Positive  ions  have  also  been  identified  (References  16  and  17)  as 
well  as  positive  cluster  ions  (Reference  18).  It  is  generally  agreed 
among  the  investigators  that  C0^  and  0^  are  the  most  dominant  positive  ions 
formed.  Since  the  positive  ions  are  not  known  to  affect  the  discharge, 
the  importance  of  negative  ions  and  negative  ion  chemistry  is  apparent. 

Exhaustive  studies  have  been  done  on  the  identification  of  con¬ 
taminants  in  a  C02/N2/He  discharge  as  evident  from  the  above.  The  only 
exception,  at  the  start  of  this  project,  was  the  identification  of  the 
negative  ion  species.  Hence,  this  identification  was  accomplished  first 
(References  14  and  15).  A  more  complete  investigation  of  the  effects 
of  hydrogen  in  a  CC^/^  discharge  was  also  accomplished  under  this 
project  (Reference  19). 

The  major  contaminant  species  CO  and  02  were  tracked  over  a  wide 
range  of  current  values  as  well  as  the  CO2  concentration.  This  experi¬ 
mental  data  was  compared  with  a  chemistry  code  (Reference  20)  in  order 
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to  identify  the  most  important  reactions  that  form  the  contaminant  species. 
The  above  comparison,  known  as  the  plasma  synthesis  approach  (Reference 

21)  was  the  last  investigation  accomplished  under  this  project  (Reference 

22) . 


The  plasma  synthesis  approach  initially  involves  the  analysis  of 
a  single  gas  component.  This  analysis  includes  a  comparison  of  the  experi¬ 
mental  results  with  the  code  calculations.  Next,  an  additional  component 
is  added  and  the  analysis  is  repeated  on  the  resultant  mixture.  This 
procedure  of  adding  components  and  then  analyzing  the  resultant  mixture 
is  repeated  until  the  total  laser  mixture  is  synthesized.  The  synthesis 
approach  is  very  easily  implemented  and  it  gives  very  reliable  results 
as  will  be  shown  later. 

The  ultimate  goal  of  this  project  was  to  investigate  the  plasma 
chemistry  of  the  C02  laser  discharge  and,  thereby,  indentifying  some 
of  the  plasma  chemical  processes  responsible  for  the  production  of 
contaminants.  A  knowledge  of  these  basic  processes  may  then  be  used 
to  make  closed-cycle  C02  lasers  more  attractive  to  the  Air  Force.  The 
Air  Force  is  interested  in  closed-cycle  lasers  not  only  for  their 
compactness,  but  also  because  they  allow  the  use  of  isotopic  C02  gas. 

The  compactness  results  from  the  elimination  of  large  storage  tanks, 
and  the  use  of  isotopic  C02  offers  the  choice  of  wavelength  selection 
corresponding  to  a  particular  atmospheric  transmission  window.  Since 
this  project  was  the  first  of  its  kind  for  this  laboratory,  experimental 
techniques  have  also  been  included,  hence,  rendering  this  report  a 
guide  in  any  future  research  in  this  area. 
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SECTION  II 

EXPERIMENTAL  ARRANGEMENT  AND  PROCEDURE 

Figure  1  is  a  schematic  diagram  of  the  experimental  arrangement. 

The  gases  used  were  of  research  grade  with  a  minimum  purity  of  99.998ft. 
Among  the  pure  gases  used  were  Ar,  He,  CO^,  N^,  and  H^.  The  gas  mixture 
used  was  1.98ft  8^/25 .3%  CO^/  72.72ft  N^.  The  mixing  manifold  (E  in  Figure 
1)  normally  contains  glass  beads  in  order  to  aid  in  a  thorough  mixing 
of  the  gas  components.  The  individual  gas  flow  rates  were  measured  by 
Hastings  mass  flowmeters.  The  pressure  transducers,  used  in  conjunction 
with  the  flowmeters,  are  represented  by  A,  B,  C  and  D  in  Figure  1.  The 
flow  rates  were  accurate  to  +  1ft  for  the  rare  gases  and  those  of  the 
molecular  gases  were  accurate  to  +  5ft. 

Several  types  of  discharge  tubes,  each  with  the  sampling  orifice 
in  a  different  location,  were  used  in  the  experiments.  The  one  shown 
as  G  in  Figure  1  has  the  sampling  orifice  (H)  on  the  side  of  the  tube 
located  in  the  middle  of  the  two  electrodes.  The  tube  was  constructed 
of  quartz  with  a  double-wall.  The  outer  jacket  was  used  for  cooling 
water.  The  inner  diameter  of  the  tube  was  2.50cm,  and  the  hollow 
electrodes  were  separated  by  a  distance  of  30cm.  The  electrodes  were 
made  of  Kovar  and  the  cathode  was  located  upstream  of  the  gas  flow. 
Pressure  measurements  were  made  behind  the  cathode  by  a  differential 
capacitance  manometer  (L)  manufactured  by  MKS  Instruments,  Inc.,  and 
the  cryosorption  pump  (M)  was  used  to  pump  out  the  reference  port  of 
the  manometer. 

The  partial  pressure  of  a  particular  gas  component  in  the  discharge 
tube  is  given  by 

» r  (3> 

where  p  is  the  partial  pressure,  p  is  the  total  pressure,  H  is  the 

5  S 

flow  rate  of  the  gas  component,  and  H  is  the  total  flow  rate.  The  total 


Figure  1.  Experimental  Arrangement  of  an  On/Off-Line  Mass 
Analyzer  System 
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flow  rate,  usually  given  in  standard  cubic  centimeters  per  minute  (SCCM), 
is  just  the  sum  of  the  flow  rates  of  the  individual  gas  components: 

H  =  E  H  (4) 

s  s 

The  gas  residence  time,  which  is  the  length  of  time  that  a  portion  of 
gas  is  in  the  discharge  before  it  is  sampled  by  the  spectrometer,  is 
given  by 


where  A  is  the  cross  sectional  area  of  the  discharge  tube,  pQ  is  standard 
atmospheric  pressure,  and  s’  is  the  length  from  the  cathode  to  the  sampling 
orifice.  Complete  derivations  of  Equations  3,  4,  and  5  are  given 
in  Appendix  A. 

A  quadrupole  mass  spectrometer,  manufactured  by  Extranuclear  Labora¬ 
tories,  Inc.,  was  used  to  sample  the  discharge.  The  spectrometer  contains 
two  differential  pumping  stages  (see  Figure  2).  The  first  stage  is  pumped 
with  a  turbomolecular  pump  and  the  second  stage  is  pumped  with  a  diffusion 
pump.  Both  pumps  are  backed  by  a  forepump.  In  normal  operation  of  the 
spectrometer,  the  pressure  in  both  chambers  are  monitored  with  ion  gauges; 
however,  during  the  pump  down  phase,  the  pressure  can  be  monitored  with 
thermocouple  gauges.  The  lowest  pressures  that  were  attainable  in  this 

experiment  were''  10"^  and"  10  ^  torr  in  the  front  and  rear  chambers, 
respectively. 

The  ion  gauge  in  the  rear  chamber  is  connected  to  an  interlock 
system  so  that  the  high  voltages  to  the  mass  filter  and  particle 
multiplier,  both  of  which  are  in  the  rear  chamber,  are  automatically 
shut  off  if  there  is  a  partial  loss  in  vacuum.  The  upper  limit  of 

the  pressure  is  pre-selected,  and  this  limit  is  usually  set  at  10'6 

torr.  For  example,  if  the  pressure  in  the  rear  chamber  exceeds  10’^ 
torr,  the  high  voltages  are  automatically  shut  off.  The  purpose  of 
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Figure  Z.  Cross  Sectional  View  of  the  Quadrupole  Mass 
Spectrometer 
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this  interlock  system  is  to  prevent  damaging  arcs  from  occurring  inside 
the  mass  filter  and  particle  multiplier  assemblies.  Because  of  the 
differential  pumping  scheme,  the  pressure  in  the  front  chamber  can  be 
-4 

as  high  as  10  torr  without  affecting  the  operation  of  the  mass  filter 
and  multiplier  in  the  rear  chamber. 

Neutral  species  were  sampled  from  the  positive  column,  while  posi¬ 
tive  and  negative  ions  were  sampled  from  the  afterglow  behind  the 
cathode  and  anode,  respectively.  The  samples  entered  the  spectrometer 
as  molecular  beams  produced  by  the  molecules  effusing  through  the  sampling 
orifice.  The  type  of  flow  from  the  orifice  can  either  be  one  of  two 
extremes:  molecular  flow  or  free  jet  expansion.  The  type  of  flow  also 
depends  on  the  size  of  the  orifice  and  the  pressure  within  the  dis¬ 
charge  tube.  The  advantages  and  disadvantages  of  each  type  of  flow 
regime  will  be  discussed  in  the  following  two  paragraphs. 

In  the  molecular  flow  regime,  molecular  species  s  flows  through 
the  orifice  at  a  rate  given  by 

<PS  =  1/4  n$  vs  A  (6) 

where  ns  is  the  steady  state  number  density  within  the  discharge  tube, 
v$  is  the  average  molecular  speed,  and  A  is  the  area  of  the  orifice. 
Molecular  flow  through  an  orifice  can  be  achieved  at  a  pressure  where 
the  mean  free  path  of  molecules  is  between  15  and  40  times  the  diameter 
of  the  orifice  (Reference  1).  The  advantages  of  molecular  flow  (also 
called  effusive  flow)  is  that  the  molecules  travel  approximately 
independent  of  each  other,  thereby  minimizing  chemical  reactions  that 
might  occur  within  the  beam.  This  aspect  of  molecular  flow  would  ensure 
that  the  chemical  constitutes  of  the  sample  are  the  same  as  those  in  the 
discharge  tube.  The  intensity  of  the  beam  falls  off  as  1 /r 2  and  is 
proportional  to  cos  (o),  where  r  is  the  distance  from  the  orifice  to  any 
point  on  the  axis  of  the  spectrometer  and  0  is  the  angle  from  the  axis 
of  the  beam  (Reference  1).  Therefore,  the  orifice  must  be  precisely 
aligned  with  the  axis  of  the  ionizer  in  order  to  obtain  a  sizeable 
signal  from  the  mass  spectrometer. 
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The  problem  of  low  sensitivity,  which  is  a  disadvantage  in  the 
molecular  flow  regime,  can  be  alleviated  by  a  free  jet  expansion  of  the 
molecules  through  the  orifice.  This  regime  occurs  when  the  mean  free 
path  is  much  less  than  the  diameter  of  the  orifice  and  mass  motion  of 
the  gas  occurs.  Also,  the  gas  in  a  free  jet  expansion  is  collision 
dominated  and  chemical  reactions  are  highly  probable.  Therefore,  the 
sample  obtained  via  a  free  jet  expansion  does  not  accurately  represent 
the  chemistry  of  the  discharge.  This  aspect  is  the  major  disadvantage 
of  a  free  jet  expansion.  A  trade-off  regime,  known  as  transition  flow, 
occurs  when  the  mean  free  path  is  on  the  order  of  the  diameter  of  the 
ori f ice. 

To  sample  neutrals,  the  molecular  beam  is  led  into  an  electron 
impact  ionizer.  The  most  common  setting  of  the  electron  energy  is 
70eV  because  the  maximum  ionization  efficiency  of  a  majority  of  the 
atomic  and  molecular  gases  occur  between  60  and  1 20eV  (Reference  23). 
Another  reason  is  that  the  energy  of  the  electrons  is  only  approximately 
mono-energetic;  consequently,  there  is  a  minimum  fluctuation  of  ion 
pair  production  from  a  fluctuation  in  the  electron  energy  in  the  region 
of  maximum  efficiency.  The  disadvantage  of  the  70eV  setting  is  that 
other  processes  such  as  double  ionization,  excitation,  and  dissociation 
have  a  high  probability  of  occurring. 

The  meter  reading  of  the  electron  energy  does  not  reflect  the  true 
energy  of  the  bombarding  electrons.  This  fact  is  due  to  the  contact 
potential  difference  of  the  electrical  connections  and  space  charge 
effects  in  the  ionizer.  The  amount  of  offset  is  usually  a  few  eV  and 
it  depends  on  the  particular  construction  of  the  ionizer.  The  offset 
can  be  determined  by  making  appearance  potential  plots  of  several  known 
gases.  The  onset  potentials  are  then  compared  with  the  published  values 
of  the  ionization  potentials.  The  difference  is  the  amount  of  offset. 

Once  the  ions  are  formed,  they  need  to  be  extracted  from  the  ionizer 
and  focused  into  the  mass  filter.  The  ions  are  extracted  by  an  extractor 
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lens  and  the  voltage  on  this  lens  can  be  varied  from  -  20  to  +  20  volts. 

To  provide  ion  focusing,  the  voltage  on  each  of  a  series  of  three  lenses 
can  be  varied  from  -400  to  +100  volts.  The  theory  of  electric  mass 
filtering  may  be  found  in  Appendix  B. 

To  detect  the  ions,  a  21 -stage  copper-beryllium  particle  multiplier 
is  used.  The  multiplier  is  positioned  off-axis  from  the  mass  filter 
in  order  to  prevent  neutrals,  photons,  and  metastables  from  striking  the 
multiplier.  Hence,  spurious  noise  is  eliminated,  and  the  signal  to  noise 
ratio  is  greatly  improved.  The  positive  ions,  therefore,  must  rely  on  the 
high  negative  potential  of  the  first  dynode  in  order  to  be  deflected 
toward  the  multiplier  and  subsequently  detected.  The  potential  of  the 
first  dynode  is  usually  kept  between  -4  and  -5kV.  When  detecting  negative 
ions,  however,  the  polarity  of  the  multiplier  must  be  reversed. 

The  signal  from  the  particle  multiplier  is  fed  into  a  preamplifier 
which  is  a  pulse  amplifier-discriminator.  Counting  or  current  measure¬ 
ments  can  be  made  simultaneously  or  individually  on  the  same  signal  source. 
The  pulse  pair  resolution  time  of  the  circuit  is  0.2  psec,  and  the  maximum 
available  count  rate  is  1  MHz.  The  count  rate  is  placed  below  the 
operating  frequency  of  the  mass  filter  so  that  no  stray  rf  can  enter  the 
circuit.  Pulse  counting  is  advantageous  because  noise  can  be  minimized 
by  adjusting  the  height  of  the  discriminator  level. 

The  signal  out  of  the  preamp  is  then  fed  into  a  fast  electrometer 
for  further  amplification.  The  detectable  current  range  of  this  instru- 
-1 1  -4 

ment  is  2  X  10  to  1  X  10  A.  The  signal  out  of  the  electrometer  is 
subsequently  displayed  by  an  oscilloscope  and/or  chart  recorder. 

It  is  highly  desirable  to  eliminate  the  background  spectrum  when 
interpreting  mass  spectra.  Modulated  beam/lock-in  amplifier  techniques 
can  be  used  to  accomplish  this  task.  A  tuning  fork  chopper  with  a 
nominal  frequency  of  400  Hz  is  used  as  the  modulator.  The  chopper 
physically  interrupts  the  molecular  beam  in  a  periodic  manner,  and  the 
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ions  arriving  at  the  particle  multiplier  have  amplitude  changes  which 
are  at  the  same  frequency  as  the  chopper.  The  ions  and  neutrals  arriving 
from  the  background  gas  have  no  such  periodic  nature.  The  lock-in 
amplifier  detects  these  amplitude  changes  and  electronically  subtracts 
the  background  signal  from  the  total  signal,  thus,  indicating  the  true 
signal . 

The  major  experimental  variables  were  the  gas  pressure,  discharge 
current,  and  gas  residence  time.  The  pressure  range  of  the  experiments 
was  from  1  to  32  torr,  and  the  current  range  was  from  5  to  80  mA.  In 
the  experiments  that  used  a  flowing  gas  mixture,  the  gas  residence  time 
was  adjusted  to  allow  many  of  the  plasma  chemical  reactions  to  reach 
equilibrium.  Equilibrium  conditions  can  be  determined  by  noting  the 
concentrations  of  chemical  species  as  a  function  of  residence  time. 

Since  the  geometry  of  the  tube  is  fixed,  the  residence  time  for  a 
particular  value  of  the  gas  pressure  can  be  varied  by  varying  the  total 
flow  rate  (refer  to  Equation  5).  Figure  3  shows  the  concentration 
versus  residence  time  of  a  few  contaminants  that  are  formed  in  a 
C02/N2/H2  discharge. 

The  concentration  in  the  discharge  can  be  predicted  with  a  high 
degree  of  accuracy  if  certain  requirements  are  met.  These  requirements 
and  the  method  for  determining  component  concentrations  are  discussed 
in  Appendix  C. 


Figure  3.  Residence  Time  Plot  (Reference  45) 
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SECTION  III 
CHEMICAL  MODELING 

As  a  first  approach  to  modeling  the  neutral  chemistry  of  the 
discharge,  the  following  simplifications  were  made: 

(1)  Nearly  all  the  reactions  considered  were  binary  and 
assumed  to  be  second  order  and 

(2)  Only  the  forward  rate  constants  were  considered,  while 
the  reverse  rate  constants  were  set  equal  to  zero. 

Considering  mostly  second  order  binary  reactions  is  reasonable  since  the 
experiments  were  carried  out  at  low  pressures  (2  TORR).  The  second 
simplification  appears  to  be  extreme;  however,  the  availability  of 
data  did  not  permit  the  inclusion  of  reverse  rate  constants.  As  will 
be  shown  in  the  next  section,  these  simplifications  (along  with  a  few 
adjustments)  resulted  in  a  chemical  code  that  gave  fair  agreement  with 
the  experimental  data. 


Consider  an  elementary  binary  reaction  with  the  following 
stoichiometry: 


a .  A.  +  b.  B .  «-  c.C.  +  d.  0. 

it  l  l  k  .  i  i  ii 

-i 


where  and  k  .  are  the  forward  and  reverse  rate  constants,  respectively. 
The  rate  of  the  reaction  is  given  by 

,  d(A  1  ,  d[B  ]  ,  d[C  )  ,  d[D  ] 

Ri  =  '  s;  -dt2-  =  '  tr  -ar-  =  ft  -dt-  =  dT  -dt-  (8) 

where  the  brackets  denote  number  density.  Since  Equation  7  is  assumed 
to  be  a  second  order  reaction,  the  rate  of  formation  of  species  Ci  is 

proportional  to  the  product  of  the  A.  and  B.  number  densities,  while 
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the  rate  of  depletion  of  is  proportional  to  the  product  of  the 

and  D..  number  densities.  The  constants  of  proportionality  are  the 

forward  and  reverse  rate  constants,  respecti vely.  Therefore,  the  rate 
of  formation  of  species  C..  is  given  by 


'c.  V1  ■  kf  [Ai'  'V  -k-t  ICi'  '“f1  -  ki  t 1  [  B,  1  (9) 


since  k  .  =  0  by  simplification  2.  can  also  be  depleted 

by  reacting  with  other  species.  For  examples,  a  depleting  reaction  is 
given  by 


ki 

a.A.+c.C.  b.B.+d.D. 

J  J  J  i  k<-  3  3  3  3 

-3 


(10) 


where  the  corresponding  rate  of  depletion  is 


.  1  dlCjJ 
c .  dt 
3 


W'i' 


-  k  . [B  .] [D .] 
-J  3  3 


=  k  .[A  .]  [C.  ] 

J  J  l 


(ID 


The  net  rate  of  formation  of  species  C..  is  determined  by  summing  over  all 
reactions  that  form  and  deplete  C.. : 


d[C]  _ 
dt 


k.  [A. ] [B. ]  -  Z 

i  1  l  1  l 1  i 
J 


c.  k.  [ A . ] [ C . ] 
J  J  J  i 


(12) 


Note  that  Equations  10  and  11  are  multiplied  by  the  stoichiometric 

factors  c.  and  c.  before  the  sum  is  taken  (which  is  the  convention 
i  3 

used  by  many  researchers).  The  result  of  this  formulation  (which  can 
easily  be  extended  to  third  order  ternary  reactions)  is  a  set  of  coupled 
differential  equations  that  are  usually  solved  numerically.  Reference 
33  gives  a  detailed  explanation  of  reactions  and  rate  laws. 
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There  are  examples  in  the  literature  of  CO^  laser  discharge  modeling. 

These  models  included  ion  reactions  as  well  as  neutral  reactions.  For 
example,  Wiegand  and  Nighan  modeled  the  CC^/N^/He  discharge  and  predicted 

that  numerous  minority  species  {such  as  those  composed  of  carbon,  oxygen, 
and  nitrogen)  are  produced  in  significant  quantities  on  a  time  scale  of 
100  ysec  (Reference  34).  These  minority  species  are  found  to  be  directly 
related  to  the  dominant  ionic  species.  Pace  and  Lacome  used  a  chemistry 
code  to  predict  that  when  small  amounts  of  0?  or  H^O  are  present  in  the 

discharge,  the  negative  ion  population  is  significantly  increased  (Refer¬ 
ence  13).  Finally,  Thoenes  and  Kurzius  have  modeled  the  plasma  chemistry 
important  to  the  formation  of  impurities  (Reference  37).  Valuable  infor¬ 
mation  concerning  the  formative  reactions  of  impurities  can  be  obtained 
from  chemical  modeling.  This  knowledge  is  valuable  because  it  often 
leads  to  solutions  for  suppressing  the  impurities. 

The  chemical  code  used  to  model  the  discharge  was  developed  by  Stamm 
(Reference  20).  The  resultant  differential  equations  were  obtained  by 
considering  the  reactions  listed  in  Table  1  and  were  solved  by  the 
modified  Euler  method.  However,  certain  assumptions  and  adjustments 
had  to  be  made  in  order  to  obtain  a  model  in  agreement  with  experiment. 
These  assumptions  and  adjustments  are  summarized  below. 

The  electron  number  density  in  the  positive  column  of  the  discharge 
is  assumed  constant,  and  it  can  be  approximately  calculated  from  the 
following  equation: 


where  I  is  the  measured  discharge  current,  A  is  the  cross  sectional  area 
of  the  discharge  tube,  e  is  the  magnitude  of  the  electronic  charge,  and 
vd  is  the  drift  velocity  of  the  electrons.  The  drift  velocity  of  the 
electrons  was  taken  from  the  data  in  Reference  38.  However,  a  knowledge 
of  the  reduced  field  in  the  positive  column  was  required  first.  It  was 
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TABLE  1 


REACTIONS  AND  RATE  CONSTANTS  USED  IN  THE  CHEMICAL  MODEL 


No. 

REACTION 

RATE  CONSTANT  (300°K) 

REFERENCE 

1. 

BINARY 

e-+C02-C0+0+e~ 

See  Text 

2. 

e_+C0-vC+0_ 

.  lr4-14  3  -1 

3x10  cm  s 

36 

3. 

e”+02->20+e~ 

1x10  ^cm^s-l 

36 

A. 

e_+N2->2N+e“ 

e~+NCHN+0" 

See  Text 

5. 

1  irT12  3  -1 

1x10  cm  s 

36 

6. 

N-KKNO 

7  1r-L7  3  -I 

2x10  cm  s 

20 

7. 

N+NO-+N2+O 

9  7  .  -11  3  -1 

2.7x10  cm  s 

39 

8. 

N+N-+N0 

q  ln-18  3  -1 

9x10  cm  s 

20 

9. 

N+0  ->N0+O 

-19  3  -] 

2.8x10  cm  s 

39 

10. 

n+no3->no+no2 

n+o3->nch-o2 

0+0-->202 

.  7  in-13  3  1 

3.7x10  cm  s 

39 

11. 

c  in-13  3  -1 

5.7x10  cm  s 

39 

12. 

c  ir.-15  3  -1 

o.5xl0  cm  s 

39 

13. 

o+no3-k>2+no2 

-13  3  -1 

5. 7x10  cm  s 

39 

1A. 

N0+N0,->2N02 

a  ,  ,  -12  3  -1 

8.7x10  cm  s 

39 

15. 

NOH)3->N02+02 

.  .  in-14  3  -1 

1.6x10  cm  s 

39 

16. 

no?+o^-+no?+o9 

3.2x10  ^cm^s  1 

39 

17. 

0+N02-K)2+N0 

no9+no^no2+o9+no 

q  7  7A-12  3  -1 

9.2x10  cm  s 

39 

18. 

.  .  in-15  3  -1 

1.1x10  cm  s 

39 

19. 

NO  3+N0  2N0  2+0  2 

2.1x10  ^cm^s  * 

39 

20. 

N+N0o-*2N0 

l.AxlO  ^cn7s  ^ 

39 

21. 

N+CO^NO+CO 

1  n  in*19  3  -1 

1.0x10  cm  s 

39 

22. 

TERNARY 

N+04-M*->N0+M 

1.0x10  ^cm^s  ^ 

39 

23. 

N+N+M-*N0+M 

,  ,  in-33  6  -1 

4  4x10  cm  s 

39 

2A. 

O+NO  +M->NO~+M 

a  £vin"32  6  “1 

9.6X10  cm  s 

39 

25. 

N+NO+M-+N  _  O-t-M 

1.0x10  ^cm^s  ^ 

20 

26. 

O+O+M-K)  +M 

7  0  in-33  6  -1 

7.2x10  cm  s 

39 

27. 

0+N  0+M-+N  0  _ +M 

q  ,  in-32  6  -1 

9.4x10  cm  s 

39 

28. 

0+0  +M->0  +M 

,  .  ln-3A  6  -1 

5. Ax 10  cm  s 

39 

29. 

o+co+M'K:o2+m 

7  .  .  -36  6  -1 

2.4x10  cm  s 

39 

*ALL 

- 

SPECIES 
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calculated  approximately  by  the  following  procedure.  The  electric  field 
was  determined  by  subtracting  the  total  tube  voltage  from  the  cathode 
fall  of  potential  and  dividing  this  voltage  difference  by  the  length 
of  the  positive  column.  The  electric  field  was  then  divided  by  the  gas 
number  density  in  order  to  arrive  at  a  value  of  the  reduced  field. 

From  the  reactions  considered,  the  rate  for  C02  dissociation  by 

electron  impact  was  the  most  dominant  term  in  the  total  rate  equation 
for  CO^.  However,  the  published  value  of  the  rate  constant  for  C02 

dissociation  (e.g.  Reference  36)  did  not  result  in  C02  dissociation  in 

agreement  with  experiment.  The  code  prediction  of  C02  dissociation  was 

far  greater  than  that  actually  observed.  A  rate  constant  providing 
agreement  with  experiment  was  obtained  in  the  following  way.  It  was 
determined  that  the  C02  rate  equation  could  be  approximated  by 

dfCO-I 

-3T-  ■  -^1  "e  <”> 

Integrating  Equation  14  and  applying  the  initial  condition  that  at 
t=0,  [CO^]  =  [C02]o,  t*ie  fol1owin9  result  is  obtained: 

-knet 

[C02]  =  [C023o  e  (15) 


Solving  for  the  rate  constant  and  using  a  time  of  1  sec  (which  is 
consistent  with  the  gas  residence  time). 


k 


i'co2]0 

[C02] 


(16) 


[C02]q  and  [C02]  can  be  obtained  from  experiment,  and  ng  can  be  calculated 

by  the  procedure  mentioned  previously.  The  rate  constant,  calculated 
by  using  Equation  16,  was  used  in  the  code  so  that  the  theoretical 
curve  for  C02  concentration  versus  discharge  current  is  made  to  fit  the 
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experimental  one.  The  concentration  of  all  other  species  were  compared 
with  the  experimental  values  with  the  above  adjustment  made  to  the  CO^ 
dissociation  rate  constant. 
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SECTION  IV 

EXPERIMENTAL  RESULTS  AND  COMPARISON  WITH  THEORY 


The  voltage-current  characteristics  of  the  mixtures  used  in  the 
synthesis  experiment  are  shown  in  Figure  4.  The  Ar  in  each  mixture 
was  used  for  calibration  purposes  and  its  concentration  was  approximately 
10%  (see  Appendix  C).  Using  the  ZO^/hr  curve  as  a  reference,  it  can  be 

seen  that  N^/Ar  and  C02/N2/Ar  discharges  operated  at  higher  sustaining 

voltages  while  CO^/He/Ar  and  C02/N2/He/Ar  discharges  operated  at  lower 

ones.  This  difference  is  explained  by  the  following  material. 


The  sustaining  voltage  of  a  glow  discharge  depends  inversely  on  the 
ionization  coefficient  which  is  defined  as  the  number  of  ion  pairs  formed 
by  an  electron  per  volt  in  the  direction  of  the  field  (Reference  40). 

(The  sustaining  voltage  also  depends  on  the  secondary  emission  coefficient 
which  appears  in  a  relatively  slow  varying  term.)  The  ionization  coefficient 
depends  on  the  nature  of  the  gas  and  the  reduced  field  of  the  discharge 
(References  23  &  40).  The  reduced  field  of  the  N2/Ar  discharge  is 

calculated  to  be  5.5  Td  using  the  method  described  in  the  previous 
section.  Consequently,  there  is  approximately  an  85%  probability 
that  a  N2  molecule  is  vibrationally  excited;  and  the  probability  of 
ionizing  a  N2  molecule  is  vanishingly  small  (see  Figure  9).  Although 
the  result  in  Figure  9  is  for  the  fractional  power  transfer  in  a  N2 
discharge,  it  can  be  translated  into  the  probability  of  occurrence 
of  the  processes  shown  as  a  function  of  the  reduced  field.  Since  N2 
constitutes  90%  of  the  mixture  and  the  ionization  of  N2  contributes 
little  to  the  ionization  coefficient,  it  can  be  seen  why  the  sustaining 
voltage  is  higher.  Also,  according  to  Figure  9,  the  formation  of  ion 
pairs  by  the  metastable  reactions  (Reference  41) 


and 


is  not  highly  probable; 


N2*  +  e"  -*■  +  2e~ 

(17) 

N2  +  e~  >  N+  +  N  +  2e~ 

(18) 

and  further  evidence  showing  that  Equation  18 
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has  a  low  probability  of  occurring  is  that  there  is  little  or  no 
dissociation  in  a  N^/Ar  discharge  (Reference  22). 

It  is  obvious  that  the  high  sustaining  voltage  of  the  CO^/N^/Ar 

discharge  is  the  result  of  a  low  ionization  coefficient.  Hence,  there 
is  very  little  ion  pair  formation  occurring.  Whether  the  discharge  is 
mostly  vibrationally  or  electronically  excited  is  unknown  because  of  the 
inavailability  of  data  such  as  that  in  Figure  9. 

When  He  is  added  to  CC^/Ar  and  CO^/N^/Ar,  there  is  a  resultant 

increase  in  the  ionization  coefficient.  This  increase  is  due  to  the 
following  Penning  reactions  (References  40  and  42): 


*  + 

He  +  Ar  +  Ar  +  e  +  He 

(19) 

He  (23s)  +  N2  -  N?+  +  e"  +He 

(20) 

As  a  result,  the  sustaining  voltage  is  lowered  considerably. 

Figure  5  shows  the  discharge  pressure  as  a  function  of  current. 

The  pressure  increased  with  current  in  all  of  the  mixtures  except  that 
of  ^/Ar.  In  the  case  of  N^/Ar,  the  pressure  dropped  by  0.01  Torr  and 

remained  constant  over  the  entire  range  of  current  values  (10-70  mA). 

This  pressure  drop  has  also  been  observed  by  Karube  and  Yamaka  in  the 
case  of  a  sealed-off  discharge  (Reference  6).  If  the  gas  mixture 
is  inert,  gas  heating  and  electrophoresis  are  the  principle  causes 
of  the  pressure  change;  and  if  the  mixture  is  chemically  active  under 
discharge  conditions,  dissociation  must  be  included  as  one  of  the  causes. 

Figures  6,  7,  and  8  summarizes  the  neutral  chemistry  of  a  CO^/Ar 

discharge.  The  average  rate  constant  (i.e.  calculated  from  Equation  16) 

11  3  1 

for  CO^  dissociation  is  5.4  X  10'  cm  s'  .  Using  this  rate  constant 
in  the  chemical  code,  the  experimental  and  theoretical  CO^  concentrations 

in  Figure  6  agree  favorably  as  expected.  The  experimental  fractional 
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dissociation  of  CO2  was  44%.  Buser  and  Sullivan  has  shown  that  dissocia¬ 
tion  up  to  80%  is  common  (Reference  5).  Dissociation  is  preceeded  by  the 
production  of  an  electronic  excited  state  of  CC^  by  electron  impact 
(Reference  7) : 

C02  +  e"  +  C02  (3*)  +  e"  (21 ) 

with  subsequent  dissociation  occurring  through  the  following  reaction: 

C02  (3tt)  +  e"  -*■  CO  +  0  +  e'  (22) 

Dissociation  can  also  occur  through  dissociative  attachment  (Reference 
36) : 


C02  +  e-  CO  +  0“  (23) 

The  overall  rate  constant  for  the  mechanism  illustrated  in  Equation  21 
and  Equation  22  is  1  X  10‘9  cm3  s"1  (3006K),  while  that  of  Equation  23 

is  5  X  10~13  cm3  s-1  (300°K)  (Reference  36).  The  calculated  rate 
constant  for  C02  dissociation  lies  between  those  of  dissociation  by 

electron  impact  and  dissociative  attachment.  This  calculated  rate 
constant  is  larger  than  that  of  dissociative  attachment  by  a  factor 

of  103. 


With  a  large  dissociation  of  C02,  there  should  be  a  large  formation 

of  CO.  This  was  indeed  the  case  as  shown  in  Figure  7.  There  should 
also  be  atomic  and  molecular  oxygen,  and  the  concentration  of  0  should 
be  larger  than  that  of  O2  since  the  rate  constant  for  CO2  dissociation 

is  much  larger  than  that  of  the  three  body  recombination  of  0  (see  Table 
1).  This  behavior  is  exactly  what  is  predicted  by  the  code  (Figure  8). 
However,  O2  was  the  only  form  of  oxygen  detected  experimentally  in  all 
of  the  mixtures. 
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The  absorption  of  atomic  oxygen  on  the  walls  of  the  discharge  tube 
is  a  possible  mechanism  by  which  it  is  lost.  Atomic  oxygen  can  also 
recombine  on  the  walls  of  the  discharge  tube  forming  0^;  hence,  the 
"effective"  rate  for  0  recombination  is  larger  than  that  shown  in  Table  1. 
This  adjustment  was  not  included  in  the  code. 

No  detectable  dissociation  of  N£  was  observed  in  a  N^/Ar  discharge. 

This  observation  is  supported  by  the  fact  that  there  was  no  N  signal 
and  no  change  in  the  signal.  The  N^/Ar  discharge  operated  at  a 

reduced  field  of  5.5  Td;  hence,  the  probability  of  vibrationally 
exciting  is  much  greater  than  that  of  dissociation  (see  Figure  9). 

The  neutral  chemistry  of  a  CC^/N^/Ar  discharge  is  illustrated  in 
Figures  10,  11,  and  12.  The  CO^  and  percentages  were  40%  and  50%, 
respectively.  The  fractional  dissociation  of  CO^  was  approximately 
63%,  and  it  was  greater  than  that  in  the  case  of  a  CO^/Ar  discharge. 

As  in  the  case  of  CO^/Ar,  there  was  no  indication  of  CO  dissociation 

at  high  discharge  currents  because  there  was  no  carbon  (C)  signal.  The 
experimental  CO  percentage  in  Figure  11  is  slightly  higher  than  that 
shown.  The  reason  is  that  CO  and  Nr,  have  approximately  the  same  mass 
number,  and  a  constant  signal  (which  was  determined  from  the  mixture 
with  the  discharge  off)  was  subtracted  from  the  total  signal  at  mass 
number  28  in  determining  the  CO  signal.  Resolution  constraints  did  not 
permit  an  accurate  determination  of  the  CO  and  signals.  However, 
there  was  N£  dissociation  (%!%),  and  it  was  deduced  from  a  rise  in 
the  atomic  nitrogen  (N)  signal  and  from  the  formation  of  nitric  oxide 
(NO).  Therefore,  subtraction  of  a  constant  signal  from  the  CO 
signal  is  not  fully  justified,  but  it  represents  a  fair  approximation. 

The  only  nitrogen  oxide  detected  was  NO,  and  the  concentration  of 
NO  stayed  at  a  constant  0.3%  of  the  total  ion  current  from  the  spectro¬ 
meter  in  the  range  from  20  to  80  mA.  Nitrogen  dioxide  (NO^)  was  not 
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Figure  9.  Fractional  Power  Transfer  in  a  N?  Discharge 


Concentration  in  a  C0o/N 


0  Concentration  in  a  C0~/N?/Ar  Discharge 


The  dotted  curve  is  the  theoretical  0  concentration. 
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formed  in  detectable  quantities,  and  nitrous  oxide  (N^O)  could  not  be 

detected  because  the  resolution  of  the  spectrometer  is  not  sufficient  in 
resolving  the  N.,0  and  CO^  signals.  The  former  conclusion  is  based  on 

the  following  fact.  A  signal  appeared  at  mass  number  46  with  no 
discharge,  and  this  signal  was  attributed  to  the  carbon  dioxide  isotopes 

017  c12  017  and  0^  C^2  0^.  If  NC^  was  formed,  then  its  signal  would 
have  increased  the  total  signal  at  mass  number  46.  However,  this  was 
not  the  case. 

In  the  case  of  the  CO^/N^/Ar  discharge,  there  is  fair  agreement 

between  the  experimental  results  and  the  results  of  the  code.  An 
average  calculated  rate  constant  of  1.7X10“^  cm2  s-^  for  CC^  dissocia¬ 
tion  was  used  in  the  model.  The  code  required  a  rate  constant  for  the 
dissociation  of  N^  by  electron  impact  (reaction  number  4  in  Table  1). 

Since  N,,  is  very  similar  in  mass  and  vibrational  energy  spacing  as  0^ 

(Reference  43),  the  rate  constant  was  assumed  to  be  the  same  as  that 
given  for  reaction  number  3.  The  code  prediction  of  the  NO  concentration 
was  less  than  5  ppm  even  at  high  discharge  currents  (%80  mA).  The 
predicted  N^O  and  NO^  concentrations  were  less  than  the  NO  concentration 
2  3 

by  factors  of  10  and  10  ,  respectively.  Since  there  was  no  NO  available 
for  calibration,  adequate  comparison  between  the  experimental  and  theoreti¬ 
cal  NO  concentrations  could  not  be  made. 

Figures  13  through  18  show  the  neutral  chemistry  of  a  (^/He/Ar 
and  CO^/N^/He/Ar  discharge.  The  fractional  dissociation  of  CO^  was 
74 %  in  (^/He/Ar  and  up  to  76%  in  (X^/N^/He/Ar .  The  corresponding 
calculated  rate  constants  for  CO2  dissociation  were  2.2  X  10’^  cm2  s"^ 
and  2.3  X  10  ^  cm2  s”\  respectively. 
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Figure  14.  CO  Concentration  in  a  CO^/He/Ar  Discharge 


EXPERIMENTAL 


theoretical  0  concentration  is  shown  as  the  dotted  curve. 
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The  theoretical  atomic  oxygen  concentration  is  shown  as  the  dotted  curve 
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In  ^/Nj/He/Ar,  small  amounts  of  carbon  (C)  were  detected,  and  the 

percentage  of  C  in  the  spectrometer  ion  current  fluctuated  between  0.09% 
and  0.26%.  The  C  may  have  come  from  the  dissociation  of  CO;  however, 
this  is  not  apparent  from  Figure  17.  The  CO  concentration,  however,  was 
strangely  larger  than  the  original  CO2  concentration  at  discharge  currents 

greater  than  20  mA.  This  is  attributed  to  an  error  -n  extrapolating  the 
spectrometer  factor  (see  Appendix  C). 

N2  dissociation  occurred  also  in  the  ^/^/He/Ar  discharge. 

Consequently,  NO  was  formed,  and  its  concentration  stayed  at  a  constant 
0.5%  of  the  total  ion  current  over  the  whole  range  of  the  discharge 
current  (10-80  mA).  There  was  no  trace  of  NO2,  and  the  presence  of  N2O 

could  not  be  detected  because  of  its  similarity  in  mass  with  that  of  CO2. 

The  maximum  NO  concentration  that  the  code  predicted  was  approximately 
2  ppm  with  N20  and  N02  at  concentrations  of  5  ppb  and  1  ppb,  respectively. 

The  ion  chemistry  and  the  effects  of  in  a  C02  laser  discharge 

has  also  been  studied.  These  results  have  been  published  elsewhere 
(References  3,  14,  15  and  19),  and  only  a  brief  summary  will  be  given 
here . 


sitive  ions  were  sampled  from  the  afterglow  behind  a  hollow 
cylindrical  cathode  in  a  12.78%  ,2Z%  He  and  12.7%  CO2/I 5.5%  N 

71.8%  He  discharge.  In  both  cases,  the  total  gas  pressure  was  2  Torr 
and  the  discharge  current  was  40  mA.  The  major  positive  ions  in  the 
Cl^/Ng/He  mixture  were  C02+,  02+,  C0+,  N2+,  C+,  N+,  0+,  and  N0+;  while 

C02+,  02+»  0+,  and  C+  were  observed  in  C^/He  (Reference  3).  Chemical 

equilibrium  of  all  of  the  positive  ion  species  was  reached  in  approximately 
10  seconds  after  the  discharge  was  turned  on. 
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Negative  ions  were  sampled  from  the  afterglow  behind  a  hollow 
cylindrical  anode  in  a  12.7%  C02/  15.5%  N2/  71.8%  He  discharge.  Changes 

in  the  relative  composition  of  the  negative  ions  were  observed  in  a 
current  range  from  0.2  to  60  mA.  The  gas  residence  time  was  62  milli¬ 
seconds,  and  the  total  gas  pressure  was  1  "orr.  The  major  negative  ions 
were  NO,,’,  N03’,  O^’,  0",  and  C03"  (Reference  15).  The  following  ions 

were  also  present:  OH’,  HC03’,  CHO",  and  trace  amounts  of  CN".  The 

dominant  negative  ion  was  N02~  while  C03’  was  observed  to  be  dominant 

only  at  low  values  of  the  discharge  current.  Chemistry  code  calculations 
indicate  C03’  as  being  dominant  regardless  of  the  current  range  (Reference 
36).  However,  if  the  nitrogen  oxides  (NO,  N^O,  and  NO^)  are  added  in  the 

code  calculations,  the  dominant  ion  would  be  N02*  in  agreement  with  what 

is  observed  experimentally  (Reference  35).  Negative  ions  can  also  be 
suppressed  by  the  addition  of  CO  to  the  discharge  (Reference  14).  The 
CO  detaches  0’  from  negative  ions  such  as  C03’,  N02’,  and  N03~. 

With  the  addition  of  2%  H2  to  a  75%  N2/15%  C02/  10%  Ar  discharge, 

the  contaminants  02  and  NO  were  reduced.  C02  dissociation  was  reduced 

by  as  much  as  20%  at  high  values  of  the  discharge  current  (80  mA)  (Refer¬ 
ence  19).  The  tradeoff  is  that  as  H2  is  added,  the  production  of  the 

contaminant  species  H2Q  increases. 
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SECTION  V 
CONCLUSIONS 

CO2  dissociation  is  the  fundamental  problem  in  closed-cycle  CO2 

lasers  because  the  dissociation  products  lead  to  the  formation  of 
contaminants.  The  dissociation  of  N2  also  aids  in  the  formation  of 

contaminants,  but  N2  dissociation  only  occurs  in  the  presence  of  CO2 

whch  acts  as  a  catalyst.  The  dissociation  energy  of  ground  state  C02 

into  ground  state  products  has  been  previously  reported  as  5.5  ev 
(Reference  5).  The  dissociation  energy  of  N2  has  been  reported  to  be 

approximately  7-10  ev  (Reference  44).  Hence,  it  is  unlikely  that  N2 
dissociation  occurs  through  the  interaction  of  N2  with  vibrationally 
excited  CO2  since  the  dissociation  energy  of  N2  is  higher  than  that 
of  COp.  The  dissociation  products  of  C02,  which  are  CO  and  0^,  are 
very  similar  in  mass  and  vibrational  energy  spacing  to  N2  (Reference  43). 
The  dissociation  energies  of  CO  and  O2  are  approximately  9-11  ev  and 
5.08  ev,  respectively  (Reference  44).  0^  would  dissociate  before  it 

acquired  enough  vibrational  energy  to  dissociate  N2  since  the  dissocia¬ 
tion  energy  of  O2  is  about  half  that  of  N2-  However,  it  is  energetically 
possible  for  vibrationally  excited  CO  to  dissociate  N2. 

C02  dissociation  occurs  by  electron  impact  and/or  dissociative 

attachment.  The  calculated  rate  constant  from  Equation  16  and  the 
experimental  data,  assuming  an  exponential  decay  of  C02,  lies  between 

the  published  values  of  the  rate  constant  for  electron  impact  dissocia¬ 
tion  and  dissociative  attachment.  In  the  four  mixtures  used  in  the 

synthesis  experiment,  the  calculated  rate  constant  differed  from  the 

-9  3  -1 

published  value  of  1  X  10  cm  s  for  electron  impact  dissociation 
by  a  factor  of  approximately  20.  However,  it  was  larger  than  the  rate 


AFWAL-TR-80-2029 


3 

constant  for  dissociative  attachment  by  a  factor  of  10  .  Therefore,  it 
can  be  concluded  that  the  dissociation  of  C02  occurs  primarily  by  electron 

impact.  Although  the  low  concentration  of  0“  would  seem  to  confirm  this 
conclusion,  it  cannot  he  used  as  exclusive  proof  since  O'  clusters  rapidly 
to  form  negative  ions  (References  14  and  15). 

The  code  predicts  that  0  is  larger  than  02  in  all  mixtures,  but  02 

was  the  only  form  of  oxygen  that  was  observed  experimentally.  The  reason 
that  the  code  predicts  such  a  high  concentration  of  0  is  because  the 
recombination  rate  of  0  to  form  02  did  not  include  wall  recombination. 

The  recombination  rate  includes  only  third  body  reaction  with  the  other 
molecules  of  the  discharge  (see  reaction  number  26  in  Table  I).  Another 
reason  why  the  code  predicts  high  0  concentrations  is  due  to  the  fact 
that  the  code  does  not  include  ion  chemistry  to  take  into  account  the 
formation  of  O'. 

In  pure  C02  or  C02/He  discharges  at  low  pressure,  the  following 
positive  ion  producing  reaction  is  evnected: 

e"  +  C02  -*  C02+  +  2e"  (24) 

The  positive  ions  0+  and  02+  are  subsequently  formed  by  fast  charge 
transfer  reactions  with  C02+  in  addition  to  electron  impact  ionization 
of  0  and  02: 


C02+  +  0  -*■  0+  +  C02 

(25) 

C02+  +  02  02+  +  C02 

(26) 

Reactions  24,  25,  and  26  suggest  that  C02+,  0+,  and  02+  would  be  among  the 
dominant  ions  formed  in  a  C02/He  discharge.  This  conclusion  is  indeed  the 
case  experimentally  (Reference  3).  In  C02/N2  and  C02/N2/He  discharges. 
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the  dominant  impurity  positive  ion  is  N0+  which  is  not  surprising  since 
NO  is  the  dominant  impurity  formed  (References  3  and  22). 

Unlike  positive  ions,  negative  ions  can  affect  the  stability  of  the 
discharge  if  their  concentration  becomes  large  (Reference  13).  CO  can  be 
added  to  the  discharge  as  a  detacher  which  limits  the  population  of  the 
negative  ions  (Reference  14).  The  trade-off  is  that  the  additional  CO 
competes  with  the  C02  for  vibrational  excitation,  and  N2  dissociation 

may  increase  by  the  interaction  of  excited  CO  with  N^.  As  a  consequence 

of  N^  dissociation,  various  nitrogen  oxides  are  produced  which  are  also 

quenchers  of  excited  C02-  Again,  since  NO  is  the  dominant  neutral 

impurity  in  a  C02/N2/He  discharge  and  0"  clusters  very  rapidly  to  NO 

forming  N02",  it  can  be  seen  why  the  latter  is  the  dominant  negative  ion. 

The  added  H2  in  a  C02  laser  mixture  intercepts  the  oxygen  produced 
by  C02  dissociation,  thereby,  limiting  the  amount  of  impurity  species 
formed.  The  most  notable  of  these  impurity  species  is  NO.  H2  also 
decreases  the  amount  of  C02  dissociation  which  is  very  desirable  for 
long-term  sealed  operation  of  the  C02  laser. 

It  is  emphasized  that  optimum  conditions  cannot  be  determined  from 
studying  the  plasma  chemistry  alone,  but  a  study  of  other  laser  para¬ 
meters  such  as  output  power  need  to  be  considered.  For  this  reason, 
there  is  a  need  for  on-line  capabilities  to  actual  laser  systems  which 
include  sampling  from  high  pressure  discharge  lasers. 
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APPENDIX  A 

RESIDENCE  TIME  AND  PARTIAL  PRESSURE  CALCULATION 

Consider  a  tube  of  circular  cross  section.  Assuming  Poiseuille  flow 
through  this  tube,  the  radial  velocity  profile  is  parabolic.  As  a  result, 
the  gas  near  the  walls  have  a  slower  velocity  and  thus  stay  in  the  dis¬ 
charge  longer.  It  is  the  gas  near  the  walls  that  the  spectrometer  samples 
because  the  sampling  orifice  is  located  at  the  wall.  In  the  following 
discussion,  it  is  assumed  that  all  the  molecules  travel  through  the  tube 
at  the  same  rate. 

The  volume  of  gas  that  passes  an  imaginary  plane  perpendicular  to  the 
flow  in  time  At  is  given  by 

AV  =  v  At  A  (A.l ) 

where  v  is  the  velocity  of  the  flow,  and  A  is  the  cross  sectional  area 
of  the  tube  which  is  assumed  constant.  Solving  for  the  flow  velocity 
in  Equation  A.l,  the  following  result  is  obtained: 

v  =  —  (A. 2) 

AtA 

AV 

is  the  volume  flow  rate  past  an  imaginary  plane  within  the  tube.  This 

flow  rate  is  usually  given  in  standard  cubic  centimeters  per  minute  (SCCM). 
A  standard  cubic  centimeter  is  the  volume  occupied  by  a  given  gas  at 
standard  temperature  and  pressure  (STP).  Standard  temperature  and 
pressure  are  defined  as  20°C  and  760  mm  of  Hg  (torr),  respectively. 

Assuming  an  ideal  gas,  the  number  density  of  molecules  in  one  standard 
cubic  centimeter  is  given  by 

nQ  =  =  2.505  X  1019  cm"3  (A. 3) 

o 

where  pQ  and  Tq  are  the  standard  values  of  pressure  and  temperature, 
respectively;  and  k  is  Boltzmann's  constant. 
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If  the  temperature  is  constant,  then  the  volume  occupied  at  any 
arbitrary  pressure  p  is  given  by 


(A. 4) 


where  AV'  is  the  volume  occupied  at  any  other  arbitrary  pressure  p'.  If 
p'  =  pQ,  then  AV'  is  the  volume  occupied  at  standard  pressure.  Substitut¬ 
ing  AV  from  Equation  A. 4  into  Equation  A. 2  with  p'  =  pQ,  the  flow  velocity 
becomes 


AV'  Po  _  H _ Po 

AtA  p  A  p 


(A. 5) 


where  H  = 


AV' 
At  - 


H  is  the  standard  flow  rate  given  in  SCCM. 


The  residence  time  is  defined  as  the  length  of  time  that  a  portion  of 
the  gas  is  in  the  discharge  before  it  is  sampled  by  the  spectrometer.  The 
residence  time  is  given  by 


T  = 


i 

V 


(A. 6) 


where  l  is  the  length  that  the  gas  travels  in  the  discharge  before  it 
reaches  the  sampling  orifice.  Substituting  the  flow  velocity  from 
Equation  A. 5  into  Equation  A. 6,  the  residence  time  becomes 

,  -  W-  (A.7) 

0 

When  the  flow  of  gas  through  the  discharge  tube  reaches  steady 
state,  the  pressure  remains  constant  and  the  flow  in  is  equal  to  the  flow 
out.  Under  this  condition  of  steady  state,  the  total  number  of  molecules 
in  the  experimental  region  of  interest  (i.e.,  from  t  point  where  it  is 
sampled  by  the  spectrometer)  is  given  by 

N  =  nQHx  (A. 8) 

where  x  is  the  residence  time  of  the  gas,  H  is  the  total  flow  rate  in 
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SCCM,  and  nQ  is  the  number  density  of  molecules  in  a  standard  cubic 
centimeter.  If  the  flow  contains  a  mixture  of  gases,  then  the  total 
flow  rate  is  just  the  sum  of  the  individual  flow  rates  of  all  the  gases 
that  constitute  the  mixture: 


H  = 


Z 

s 


H 


s 


(A. 9) 


The  total  number  density  is 


n  =  N  =  n0  h  t 

Ait  A  «, 


(A. 10) 


Similarly,  the  number  density  of  species  s  is  given  by 


n  H 

n  =  --9-  -!■  T 
s  A  l 


(A. 11) 


The  fractional  amount  of  species  s  is  obtained  by  dividing  A. 11  by 
A. 10: 


Finally,  assuming  an  ideal  gas,  the  partial  pressure  of  species  and 
the  total  pressure  is  proportional  to  the  number  density  of  species  s 
and  the  total  number  density,  respectively: 


P 

P 


s 


(A. 13) 
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APPENDIX  B 

THE  ELECTRIC  MASS  FILTER 

Paul  and  Steinwedel  (Reference  24)  in  Germany,  and  R.F.  Prost  in  the 
United  States  proposed  that  ions  of  different  charge  to  mass  ratios  can 
be  separated  by  a  high-frequency  quadrupole  electric  field.  This  method 
was  tested  (Reference  25),  and  a  more  accurate  investigation  was  done 
(Reference  26).  The  following  discussion  represents  a  summary  of  the 
important  aspects  of  a  mass  filter. 

The  mass  filter  is  constructed  ideally  of  four  parallel  electrodes 
of  hyperbolic  cross  section  in  a  square  array  (refer  to  the  top  of  Figure 
19).  If  a  potential  +  is  applied  to  one  set  of  opposite  electrodes 

and  -  4>  to  the  other  set,  a  hyperbolic  electric  potential  between  the 

electrodes  is  the  result.  The  potential  is  symmetric,  two  dimensional 
(i.e.  independent  of  z),  and  the  potential  on  the  z  axis  is  zero.  The 
potential  is  described  in  general  by 

4>  =  4>q  (a  x2  +  By2)  (B.l ) 


where  a  and  B  are  constants,  and  4>o  is  the  magnitude  of  the  potential 

given  to  the  electrodes.  4>o  =  U  +  V  cos  w  t  where  U  is  a  D.C.  voltage, 

V  is  the  peak  amplitude  of  a  radio  -  frequency  (rf)  voltage,  and  oj  =  2lif 
where  f  is  the  frequency  of  the  rf  voltage.  In  practice,  the  potential 
is  approximated  by  replacing  the  hyperbolic  electrodes  with  cylindrical 
electrodes  (shown  at  the  bottom  of  Figure  4).  This  approximation  is  good 
in  the  region  near  the  axis  of  the  mass  filter. 


2  2 

In  the  space  between  the  electrodes,  V  4>  =  0  so  that  V  $  =  24>o  (a  +  B) 
=  0.  Since  4>Q  f  0,  then  B  =  -a.  The  two  boundary  conditions  are  given 
by 

4>  =  4>0  at  (xQ,  0)  and  (-xQ,  0) 


8 


IDEAL  CONSTRUCTION  WITH  ELECTRODES  OF 
HYPERBOLIC  CROSS  SECTIONS 


PRACTICAL  CONSTRUCTION  WITH  ELECTRODES  OF 
CIRCULAR  CROSS  SECTIONS 


Figure  19.  Electrode  Configuration  of  the  Electric  Mass  Filter 
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4>  =  -4>0  at  (0,  yQ)  and  (0,  -yQ) 


where  x  and  y  are  the  distances  from  the  z-axis  to  the  electrodes, 
o  o 


Since  the  electrodes 


Applying  the  two  boundary  conditions,  a  =  \  and  \ 

xo  ^o 

are  in  a  square  array,  x  =  y  =  r  where  r  is  the  radius  of  the  circle 

OOO  0 

centered  on  the  z-axis  that  is  tangent  to  the  surface  of  the  electrodes. 
The  final  form  of  the  potential  is  given  by 


2  2 

$  =  (U  +  V  cos  cut )  - — 

ro 


(B.2) 


From  Equation  B.2,  the  equations  of  motion  for  a  single  charged  ion  are 


m  z  =  0 


tot )  —7-  -  0 

(B.3) 

r0 

cot)  =  0 

(B.4) 

ro 

(B.5) 

The  D.C.  field  has  a  focusing  effect  on  the  ions  in  the  x-  direction 
since  x  *  -  Ux  and  a  defocusing  effect  in  the  y-  direction  since  y  Uy. 

The  z-component  of  the  velocity  is  constant. 

Defining  <ot  =  2  a  Equations  B.3  and  B.4  are  transformed  into 


d2x 

— j  +  (a  +2q  cos  2p)  x  =  0 
dp^ 


(B.6) 


-  (a  +  2q  cos  2p)  y  =  0 
dp 


(B.7) 


8eU  4eV 

where  a  =  — ^  2  anc*  q  =  — 2  2'  Equations  B.6  and  B.7  are  in  the  form 

mr  to  mr  (1) 

0  0 
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of  Mathieu's  differential  equation.  The  solutions  are  stable  if  they 
tend  to  zero  or  remain  bounded  as  p  approaches  +  00  and  unstable  if  they 
tend  to  +  «  as  p  approaches  +  <*>.  The  stable  solutions  are  known  as  the 
Mathieu  and  modified  Mathleu  functions,  respectively;  and  the  path  of  an 
ion  is  a  high  frequency  oscillation  about  the  z  -  axis.  Ions  that  have 
stable  trajectories  are  able  to  pass  through  the  filter  to  the  detector. 

The  unstable  solutions  are  exponential  in  nature,  and  the  path  of  the  ion 
diverges  exponentially  away  from  the  z-axis  until  it  eventually  strikes 
an  electrode  and  become  neutralized.  Stability  depends  on  the  choice  of 
values  for  the  dimensionless  parameters  a  and  q.  For  a  detailed  explana¬ 
tion  of  Mathieu  functions  and  stability,  see  Reference  27. 

Figure  20  shows  the  range  of  q-values  corresponding  to  a  certain 
a-value  for  which  the  solution  is  stable.  An  operating  line  is  established 
by  a  constant  a/q  ratio,  and  this  can  be  set  by  choosing  the  ratio  U/V. 

If  the  operating  line  intersects  the  stability  region  (e.g.,  at  points 
A  and  B  in  Figure  20),  then  there  are  (a,q)  coordinate  pairs  on  this 
line  in  the  stability  region  that  will  produce  a  stable  solution. 
Consequently,  ions  belonging  to  a  certain  e/m  range  have  stable 
oscillations  since  a  and  q  depend  on  e/m.  As  the  operating  line 
approaches  the  apex  (point  c)  of  the  stability  region,  the  e/m  range 
becomes  smaller,  and  the  number  of  ions  that  are  able  to  pass  through 
the  filter  becomes  less.  When  the  operating  line  touches  the  apex, 
the  e/m  range  is  reduced  to  zero  (infinite  resolution)  and  there  is  a 
no  transmission  of  ions. 

The  a  and  q  values  for  infinite  resolution  are  0.237  and  0.706, 
respectively.  Using  the  definitions  of  a  and  q,  and  assuming  ions 
with  a  single  charge. 


U  =  1.212  M  f2  r20 

(B.8) 

V  =  7.219  M  f2  r2 

0 

(B.9) 

where  M  is  given  in  a  m  u,  f  in  MHz,  and  rQ  in  cm.  The  numerical 
constants  in  Equations  B.8  and  B.9  are  not  unique,  but  depend  on  the 
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values  of  a  and  q.  The  ratio  U/V  has  a  value  of  0.16784  in  the  infinite 
resolution  mode.  Under  normal  operating  conditions,  therefore,  the 
resolution  can  be  adjusted  by  varying  the  two  voltages,  but  the  ratio 
U/V  must  be  kept  below  0.16784  for  adequate  transmission  of  ions.  For 
high  resolution,  U  and  V  must  be  chosen  so  that  the  operating  line 
intersects  the  stability  region  near  the  apex. 

For  mass  separation  to  occur,  the  ions  of  the  incorrect  e/m  must 
remain  in  the  filter  long  enough  so  that  they  will  be  rejected.  Because 
of  this  requirement,  the  electrodes  must  be  a  minimum  length  for  ions 
of  a  certain  energy.  If  the  ions  are  injected  into  the  filter  too  far 
off  from  the  z-axis,  they  will  strike  the  electrodes  even  if  they  have 
the  correct  e/m.  Hence,  there  is  a  restriction  on  the  maximum  distance 
from  the  z-axis  that  ions  can  be  injected.  The  filter  length  and  injec¬ 
tion  radius  of  the  ions  influence  both  the  resolution  and  transmission 
(Reference  26). 

From  Equations  8.8  and  8.9,  mass  scanning  can  be  accomplished 
by  either  varying  U  and  V,  or  by  varying  f.  If  U  and  V  are  varied, 
they  must  be  varied  by  the  same  amount  in  order  to  maintain  a  constant 
U/V  ratio  and,  hence,  constant  resolution.  All  practical  quadrupole 
power  systems  are  designed  so  that  U  and  V  are  varied  while  f  is  held 
constant.  Also  from  Equations  B.8  and  B.9,  the  maximum  mass  range 
which  can  be  scanned  depends  on  f,  U,  V,  and  rQ.  rQ  can  be  varied  by 

varying  the  diameter  of  the  electrodes.  The  maximum  values  of  U  and  V 
depend  on  the  high  power  amplifier,  and  the  maximum  mass  range  depends 
only  linearly  with  U  and  V.  However,  the  maximum  mass  range  depends 
on  the  square  of  f  and  rQ  so  that  a  large  change  in  mass  range  requires 

only  a  minor  change  in  the  system's  parameters  f  and  r  . 
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APPENDIX  C 
DATA  REDUCTION 

In  order  to  arrive  at  a  prediction  of  component  concentrations, 
an  expression  for  the  total  ion  current  density  formed  by  the  ionizer 
must  be  derived  first.  Figure  21  shows  a  schematic  diagram  of  the 
axial  ionizer  most  commonly  used  in  an  Extranuclear  model  EMBA  II  mass 
spectrometer.  The  neutral  beam  is  admitted  into  a  cylindrical  grid 
which  is  surrounded  by  four  filaments  arranged  in  a  square.  Once 
inside,  ions  are  produced  by  the  interaction  of  neutral  molecules 
with  electrons  of  a  certain  energy  formed  by  emission  from  the 
filaments  and  subsequent  acceleration  by  the  positive  voltage  on  the 
grid. 


The  number  of  ionizing  collisions  per  unit  volume  per  unit  time 
between  electrons  and  molecules  is  given  by  (Reference  28): 

n  =  ne  \  <  Vas  (V)  >  (C.l ) 

where  ng  =  the  electron  number  density, 

n^  =  the  molecular  number  density  of  species  s, 

V  =  the  magnitude  of  the  relative  velocity, 

os  (V)  =  the  total  ionization  cross-section  of 
molecular  species  s,  which  is  a 
function  of  the  relative  velocity. 


The  <  >  denotes  a  calculated  average  in  velocity  space  of  the 

product  Vos  (V).  The  electrons  possess  most  of  the  kinetic  energy 

so  that  V  2  v  e  where  v  is  the  speed  of  the  electrons.  This 

conclusion  Is  justified  as  follows. 
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Figure  21.  Axial  Electron  Impact  Ionizer 


For  molecular  flow,  the  average  speed  at  which  the  molecules  travel 
through  the  orifice  is  given  by 


v 


s 


8kT  Z2 
Trm 

s 


(C  .2) 


where  T  is  the  temperature  of  the  source,  m$  is  the  molecular  mass,  and 

k  is  Boltzmann's  constant.  This  expression  for  the  average  speed  assumes 
that  the  velocity  distribution  of  the  molecules  in  the  source  is  Maxwellian. 
For  a  source  containing  Ar  at  300° K,  v$  ~  398  m/sec.  In  a  free  jet 

expansion,  the  terminal  speed  of  the  molecules  when  they  reach  the 
ionizer  is  given  by 

*s  ■  I  (  )  <  51  )  l15  (C.3) 

y  -1  s 

where  y  is  the  ratio  of  the  specific  heat  at  constant  pressure  to  the 
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specific  heat  at  constant  volume  (Reference  29),  and  for  an  ideal  monatomic 
g as  y  =  5/3.  Using  the  same  Ar  example,  v$  ~  558  m/sec.  From  the  examples 

above,  the  range  of  speeds  are  between  398  and  558  m/sec  corresponding  to 
the  three  types  of  flows:  molecular,  transition,  and  free  jet  expansion. 
The  most  common  setting  of  the  electron  energy  is  70eV  for  reasons 
explained  in  Section  II.  The  speed  of  the  electrons  corresponding  to 
this  value  of  the  energy  is  approximately  4.96  X  10^  m/s ec.  Thus,  the 
electron  speed  is  approximately  890  times  the  speed  of  the  molecules 
in  a  free  jet  expansion. 

Replacing  the  magnitude  of  the  relative  velocity  V  with  the  speed 
of  the  electrons  vg,  Equation  C.l  reduces  to 

h  =  n  v  n'0.  (C.4) 

e  e  b  s 

The  number  of  ions  formed  per  second  in  the  ionizer  is  obtained  by 
multiplying  Equation  C.4  by  the  volume  of  the  interaction  region. 

Hence, 

N  =  Ain  v  on'  (C.5) 

e  e  s  s 

where  A  is  the  cross  sectional  area  of  the  entrance  aperture  to  the 
grid  and  L  is  the  length  of  the  grid  (see  Figure  21).  If  these 
ions  are  now  swept  out  as  fast  as  they  are  produced,  then  the  ion 
current  density  becomes 


=  en  v  La  n' 

6  6  S  S 


(C.6) 


In  Equation  C.6,  the  assumption  has  been  made  that  the  molecules 

are  singly  ionized.  Also,  the  assumption  has  been  made  that  the 

Ionization  process  is  not  having  to  compete  with  other  processes 

such  as  inelastic  scattering,  dissociation,  etc.  The  product  en  v 

6  6 
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is  the  electron  current  density;  thus,  Equation  C.6  can  further 
be  reduced  to 


(C.7) 


For  a  detection  system  that  has  a  linear  response,  the  signal 
is  proportional  to  the  ion  current  density.  Hence, 


(PH)S 


C  j  La  n 
s  e  s  s 


(C.8) 


where  ( PH ) s  is  the  "peak  height"  from  the  recorded  mass  spectrum  and 

Cs  is  a  constant  of  proportionality.  The  constant  of  proportionality 

is  in  general  different  for  different  species  because  they  do  not 
all  give  the  same  response  when  striking  the  multiplier. 


The  number  density  in  the  ionizer  in  terms  of  the  number  density 
in  the  source  in  the  case  of  molecular  flow,  is  deduced  as 


n  = 
s 


1 


4trD 


A  n 
o  s 


(C.9) 


where  n  is  the  number  density  in  the  source,  A  is  the  cross  sectional 
s  o 

area  of  the  sampling  orifice,  and  D  is  the  distance  from  the  orifice 

to  the  ionizer  (Reference  1).  Substituting  Equation  C.9  into  Equation 

C.8,  the  following  result  is  obtained: 

(PH)S  =  (CsjeLAo/4TTD2)  0sn$  (C.10) 

2 

The  spectrometer  factor,  Fs>  is  defined  as  CsjeLAo/4iTD  .  Therefore, 

{PH)s  =  Fsasns  (CJ1) 


The  spectrometer  factor  depends  on  system  parameters.  It  depends 
on  the  characteristics  of  the  detection  system  through  the  constant 
of  proportionality  C$.  It  depends  on  the  characteristics  and  geometry 
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of  the  ionizer  through  the  factors  j  and  L,  respectively.  It  also 

depends  on  the  geometry  of  the  sampling  orifice  through  the  quantity 
Aq  and  the  geometry  of  the  mass  spectrometer  through  the  quantity  D. 


The  units  of  Fg  depend  on  the  units  used  for  the  peak  height.  If 
the  peak  height  is  read  from  a  previously  recorded  mass  spectrum  the 

? 

units  are  usually  given  in  centimeters.  Hence,  the  units  of  F$  are  cm  . 

On  the  other  hand,  if  the  peak  height  is  recorded  on  a  voltmeter  connected 

to  the  electrometer  or  lock-in  amplifier,  the  units  are  given  in  volts. 

Therefore,  the  units  of  F  in  this  case  are  Vcm. 

s 

There  are  three  primary  requirements  of  gas  analysis  by  mass 
spectrometry:  (1)  the  peak  height  should  be  directly  proportional 
to  the  number  density  of  the  gas  sample  in  the  source,  (2)  in  a 
gas  mixture,  the  sum  of  the  pe'-  heights  should  be  a  linear  combination 
of  the  individual  peak  heights,  and  (3)  the  gas  flow  into  the  spectro¬ 
meter  should  remain  constant  during  a  run.  All  of  these  requirements 
are  satisfied  in  the  molecular  flow  regime  (Reference  1). 


Comparing  requirement  (1)  with  Equation  C.ll,  the  product  F$os 

is  a  constant  for  any  particular  component.  Figure  22  shows  a  plot  of 
the  peak  height  versus  number  density  in  the  case  of  flowing  the  atomic 
gas  Ar  and  the  molecular  gas  individually  through  the  discharge  tube. 
The  linearity  of  these  curves  are  clearly  evident,  and  the  slopes  of  the 

a. 


lines  are  equal  to  F^cj^  and  F^ 


^  for  Ar  and  Nj,  respectively. 


For  a  gas  mixture,  the  spectrometer  factor  of  a  particular  component 
will  differ  from  the  value  it  has  in  the  case  of  individual  gas  flow 
through  the  discharge  tube.  This  difference  can  be  explained  by  con¬ 
sidering  requirement  (2).  In  mathematical  terms,  requirement  (2)  can 
be  expressed  as 


2  (PH) 
s  s 


Ik  F  a  n 
s  s  s  s  s 


(C.12) 
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The  signal  contribution  of  a  particular  gas  component  is  determined 
by  tuning  the  mass  filter  to  the  mass  number  of  the  component.  Equation 
C.12  then  reduces  to 


(PH)S  -  ksFsosns  =  F$'osns 


(C.13) 


where  the  coefficient  k$  has 


been  included  in  the  spectrometer  factor 


Figure  23  shows  a  plot  of  the  peak  height  versus  number  density 
for  Ar  and  N ?  in  the  case  of  an  Ar/N^/CC^/He  mixture  flowing  through 

16  -  3 

the  discharge  tube.  Note  the  linearity  below  15  X  10  cm  ;  however, 

16  3 

the  slopes  are  different  from  the  ones  in  Figure  22.  At  27  X  10  cm’  , 
the  deviation  of  the  points  from  linearity  is  due  to  the  fact  that 

18  _  3 

the  total  number  density  in  the  tube  exceeded  10  cm  and  molecular 
flow  was  no  longer  valid  for  the  si2e  of  orifice  used.  The  diameter 
of  the  orifice  was  2  mils  in  both  individual  and  gas  mixture  flow. 

In  general,  the  spectrometer  factors  of  the  components  in  a 
gas  mixture  should  be  determined  experimentally  by  using  Equation  C.13 
and  known  number  densities.  This  step  must  be  done  before  the  dis¬ 
charge  is  turned  on.  The  factors  are  then  used  for  normalizing 
component  concentrations  in  the  discharge.  For  new  species  produced 
in  the  discharge,  an  identification  is  made  first;  and  then  the  spect¬ 
rometer  factors  of  these  new  components  are  determined  by  the  method 
mentioned  previously. 

Under  discharge  conditions,  there  will  be  temperature  and  density 
variations  in  the  tube  due  to  gas  heating.  Therefore,  the  signal 
must  be  corrected  in  order  to  reflect  the  true  signal  (hence,  true 
number  density)  in  the  absence  of  the  effects  of  gas  heating.  There 
is  a  method  for  calculating  the  temperature  rise  and  radial  profiles 
in  typical  C02  laser  mixtures  (Reference  30).  This  method  is  somewhat 

laborious  and  requires  the  knowledge  of  specific  parameters  concerning 
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(UJO)  1HDI3H  )Hb3d 


Figure  23.  Signal  Response  in  the  Case  of  Gas  Mixture  Flow 
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the  discharge.  To  monitor  density  changes  due  to  gas  heating,  Ar  is 
added  to  the  mixture  in  question  with  a  concentration  less  than  or 
equal  to  10%.  Since  Ar  is  inert,  the  change  in  its  peak  height  is 
due  to  density  variations  alone.  Hence,  the  peak  height  is  corrected 
by  multi  pi  ing  it  by  K^,  where  K^r  is  defined  as  the  ratio  of  the  peak 

height  with  the  discharge  off  to  that  of  the  peak  height  with  the 
discharge  on.  K^r  is  also  used  to  correct  the  other  signals  for 

the  effects  of  gas  heating.  The  assumption  is  that  all  gas  components 
are  heated  uniformly  throughout  the  volume.  This  assumption  is  not 
rigorously  correct,  but  it  represents  a  fair  approximation. 

It  must  be  emphasized  that  the  argon  correction  factor  attempts 
to  correct  for  the  effects  of  gas  heating  alone.  Pressure  changes 
when  the  discharge  is  on  may  also  be  due  to  electrophoresis  (References 
31  and  32)  as  well  as  chemical  dissociation. 
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